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Abstract 



Propagation of the heavy quark in the field of a static antiquark source is studied 
in the framework of effective Dirac equation. The model of QCD vacuum is described 
. by bilocal gluonic correlators. In the heavy quark limit the effective interaction 

\ is reduced to the potential one with 5/6 Lorentz scalar and 1/6 Lorentz vector 

^j" ' linear confinement, while spin-orbit term is in agreement with Eichten-Feinberg- 

■ Gromes results. New spin-independent corrections to the leading confining regime 

\ are identified, which arise due to the nonlocality of the interaction in time direction 

<3\ \ and quark Zitterbewegung. 

Qh! a lot of evidence exists that the non-abelian nature of Yang-Mills QCD leads to the 

Q^l confinement of colour charges. Apart from purely theoretical considerations, the main 

bulk of data on confinement comes from the lattice QCD simulations and from the phe- 
nomenology of hadronic spectra. The lattice calculations firmly establish the linear rising 
. force between two static colour sources, while the hadronic masses are most successfully 

\ described with the effective qq potential which is the sum of linear and Coulomb forces. 

a \ Complementary to these facts is the idea that QCD at large distances is a string the- 

ory, and linear potential between heavy constituents is a manifestation of the string-type 
dynamics. 

In general, the dynamics governed by QCD should be nonlocal; nevertheless, it is 
natural to assume that in the heavy quark limit it is reduced to the nonrelativistic local 
potential acting between quark and antiquark which is supplied by subleading 0{l/m?) 
corrections. Among these corrections an important role is played by spin-dependent 
forces which define the spin splittings of heavy quarkonia and serve as a testing ground 
for various theoretical approaches. 
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The most consistent derivation of the spin-dependent potentials was performed in the 
framework of the Wilson loop approach |T|, where the potentials were expressed in 
terms of expectation values of gluonic field insertions into a Wilson loop corresponding to 
the propagation of qq pair, and fundamental relations between static and spin-dependent 
potentials were established 0. 

One of the challenging problems is the long range part of the spin-orbit force 

Vsoir) = f ^ - ^] (-— + ^— (a I - (1) 

which is sensitive to the Lorentz nature of the static confining interaction e{r). If it is a 
local Lorentz scalar potential, one has 

V, = -e V2 = 0, (2) 

while for the time component of the Lorentz vector it is 

V^ = V2 = e. (3) 

Phenomenological analysis of heavy quarkonia spectra clearly prefers possibility 
but it was shown within Vacuum Background Correlators Method [3-6] and in the frame- 
work of the Coulomb gauge QCD Hamiltonian that relation (H) is respected at large 
distances without ad hoc assumption of scalar confinement. It was demonstrated in [3-6] 
that nonrelativistic reduction of the long range interaction is not the whole story, and 
additional contributions to the spin-orbit force exist, which are due to the nonlocality of 
the QCD generated interaction. 

The aim of the present paper is to study systematically the nonlocality corrections 
to the effective long range interaction for heavy quark. We have found a new constraint 
on the parameters of the interaction, which appears to be crucial for the selfconsistent 
potential-type dynamics of heavy quarks. We adopt the approach suggested recently in 

and, in a more simple version, in [0, that allows, as a byproduct, to establish explicitly 
the Lorentz nature of confinement. Our derivation is restricted to the case of a heavy 
quark propagating in the field of an infinitely heavy antiquark colour source, and is only 
the first step towards defining the full dynamics of heavy quarkonia. 

The starting point of approach [^ is the Green function SqQ for the qQ system, written 
in the Euclidean space as 

SgQ{x,y) = ^J DijDij+DA^exp |-i J d^xF^l - J d'x^+{-zd - im - x (4) 

xi)^{x)SQ{x,y)i){y), 

where Sqi^x^y) is the propagator of the static antiquark placed at the origin. To consider 
the one-body limit it is convenient to choose the modified Fock-Schwinger gauge ||10[ 



^4(^4, 0) = 0, xA{x4:, x) = 0, (5) 
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in which S'q(x, y) is simply 

Sq{x, y) = i^^d{x, - y,)e-^'^^*-y*^ + i^^e{y, - x,)e-^'^y--^-\ (6) 
Integration over gluonic field in can be performed with the result 

SqQ{x,y) = DipDip^ew\^- J d^xLeff{iJ,iJ^)^iJ^{x)SQ{x,y)ij{y), (7) 

where Leff{'ip,ip~^) is the effective quark Lagrangian: 

J d^xL,fj{ilj,^^) = J d^x^P+ix)i-td-im)^''{x) + J d^x^P^ixh^^lj^ix) < A^)^ > + 

+ d'x,d^X2^+^{xih,y'{x,)^^^{x2h,,^^'{x2) < A^.^i^M.^^M) >+•••' (8) 

where all a-s and (3-s are fundamental colour indeces, and the irreducable correla- 
tors < 74^^^^^(xi) . . . y4^„^"(x„) > of all orders should enter. The first one, < A^'^ >, 
vanishes due to the gauge and Lorentz invariances, and in what follows we keep only bilo- 
cal correlator < Afj_'^{x)A^J{y) >= JC^uJsi^yV) disregard the contributions of higher 
correlators. 

Using gauge invariance of the vacuum one has 

/C,.^J(^,1/) = (AJ^(A,)]/C^t(^,2/) = 2{XXp{X,)jK,^{x,y), (9) 

and, because of the relation {Xa)'^{Xa)] = ^^s^p ~ W^^P^s^ expression (|D takes the form 
J d'^xLeff{i).,ip^) = J d'^x4)^{x){—id — im)^°'{x) + 

+ d^xd'yi^t{x)l,i^^{x)i^;{y)lur{y)K,.{x,y) (10) 
in the limit Nq oo, yielding the Schwinger-Dyson equation 

- im)S{x, y)-ij d'^zM{x, z)S{z, y) = d'-^^ {x - y) (11) 
with the self-energy part M(x, z) given by 

- iM{x, z) = Kf,^{x, z)j^S{x, z)-fu, (12) 

where S{x, n) = ^ < i^'^ {x)ip^ (y) > is the colour trace of the quark Green function. As 
in gauge (^ Green function (|D of the static source is unity in the colour space, quantity 
S{x,y) completely defines propagation of the colourless qQ object. 
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Gauge condition (|]) can be rewritten as 

A2{x4,x)= / F^^{x4, ax)da, (13) 



A't{xi,x)= [ axkFI^^{x4,,ax)da, i = 1,2,3, (14) 
Jo 

so the average K^^, can be expressed in terms of field strength correlator < F^ui.^)F^p{y) >, 
for which we use the parametrization ^: 

where the second term A^^^^ is a full derivative and does not contribute to the confinement. 
As we are interested only in long range force, we consider only the term proportional to 
D{x — y) in (p!5|), which, in contrast to A^^), contributes to the area law with the string 
tension 

a = 2 dr dXD{T,X). (16) 
JO Jo 

Function D{u4, \u\) is actually a function of u1 + due to Lorentz invariance, but in 
our apparently non-invariant treatment we keep dependences on \u\ and 1^4 1 separately 
as in ([T6|). 

Finally, for average K^^{x,y) one has {t = x^ — 7/4): 

Kuir, X, y) = (xy) Jq da dpD{T, \ax- (3y\), 

Kii{T, X, y) = Kii{T, X, y) = 0, (17) 

Kikir, X, y) = {{xy)5ik - yiXk) ada j3d(3D{T, \ax- j3y\). 
The system of equations ([TT|), ([T^) may be rewritten in terms of wave functions as 

(-i4 - im)ip{x) + / d'^zKi^„{x, z)-f^S{x, z)-f^ip{z) = 0. (18) 



This equation is essentially nonlinear, since the eigenf unctions ipn enter the spectral 
representation for S{x, z). In the heavy quark limit we solve it perturbatively, substituting 
the free Green function Sq{x,z) into the self-energy part. The resulting linear equation 
was considered in 1|, |^ . 

As correlators (p!5| ) are defined in the Euclidean space, it is convenient to formulate 
the eigenvalues problem for equation (|l^) in the Euclidean space too and perform the 
Wick rotation to the Minkowski one afterwards arriving at the Dirac-type equation 

(a^+7om + 7oM)V' = Et/-, (19) 



4 



with operator M given by 



M{x, z) = -i drK^^ir, x, z) \ yr^e^'^^'-^^x 



/o 



(2vr) 



X 



where e = £(A;) = yk'^ + m? and 7-matrices are Euclidean ones (74^ = 70M, 7b = —ilM)- 
Operator M is nonlocal both in space and time, and our strategy is to find the leading 
local limit, to establish the first order nonlocal corrections and to check if they are small. 
To this end we approximately rewrite the expression in the curly brackets in ( pOD in the 
form 




introducing the quarks binding energy Eq = E — m. 

The leading local confining interaction is obtained after omitting the terms propor- 
tional to ik and k"^ in (0), and yields 

Mo{x, z) = 5(3) - ^) drK^^ir, f , x) {7m^^7. + l^.^^luC-'"'^ 

= 6^^\x-z)V,onf{x). (22) 

It is relevant now to comment on the Lorentz nature of confinement and to clarify 
some related confusions. The underlying interaction is bilinear in vector verteces, as it is 
clearly seen from the effective Lagrangian (H): everywhere including resulting expression 
the interaction contains 7^ ... 7,^ product. As it will be shown below, for heavy quarks 
this structure actually reduces to 70 ... 70 product, and we agree at this point with the 
statement made in the interaction is time-like vector one. However, this almost trivial 
observation does not straightforwardly help to answer another question, which is usually 
put in connection with the problem of the Lorentz structure: the effective local confining 
interaction Mq{x, z) can be proportional either to unity or to 70, and it is added either 
to the mass term (scalar confinement) or to the energy term (vector confinement) in the 
effective Dirac equation ([T^). The answer to this question depends on the behaviour of 
the function K^y. 
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The correlator D{u4, \u\) should decrease in all directions in the Euclidean space so 
that the string tension (|TBp were finite, and the Dirac structure of the confining potential 
depends on the correlation length Tg which governs this decrease (r = |x|): 

K„..f(r) = r(l + 74) / dT I d\D(r, A)U--+I?-- + ^I + 



3 r ?>r'^ 



+ r(l - 74) /„ dr I dXDir, A) - ^ + ^ + (l - ^) e'^^^ (23) 



with large distance (r ^ Tg) behaviour 

W)=(^i7.).. + 0(^) (24) 

for mTg ^ 1, and 

5 

Vconfir) = -ar + O (armTg) (25) 

for mTg <C 1. Both regimes were established in @, 0. In what follows we demonstrate 
that only the regime mTg ^ 1 is self consistent. 

Now we consider various corrections to the leading interaction (^). The terms pro- 
portional to (7^) in ( pT|) are calculated by means of integration by parts and give 

M.^.ix, z) = 5^'\x - z)V.^,{x), (26) 

i { (9 1 

%(^) = ~ 2m io ^'^^^^ e"^™^) I 7/.7i7:.^Mi^(^, ^, x)pi - i'Ji.'ya^-g^K^^iT, f, z)\g=£ \ = 

1 1"'' ( -~ / 2 X'^ \ 

= -J^ dT{l + e''-nl rfAD(r, A) I (r - A)2(7P) + f-r- A + —jz(7n)(np) + 

^ X d 
n=-, p=-t— 
r ox 



Problems start with the terms proportional to k"^ in ([211). To bring these terms into 
local form one should integrate by parts twice: 

M,,(f, z) = 6'~'\x - z}V.,{x), (27) 



2 



1 ( B 

Vp{x) = dT^,{l + + (1 + 74)mr)7. q^K^^t, x, z}y=,+ 
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d -^2] 
+2i—K^^{t, X, z)\^=xpi - K^^{t, X, x)p ^ , 

and only momentum dependent terms can be expressed in terms of integrals of the function 
D(r,A): 



1 



Jo 



dT / dXD{T, A) (1 + e"^"" + (1 + 74)mr 



+ (1 + + (1 - 74)mr) 

_^/2_A A^^ 
2r V 3 r 3r'^ 



1 /2 



A3 \ 



1 - i - 

-(r - A)p^ + -(np) 



A3 



2 \ 3 



[al) 



3r2 



+ V"^2 (momentum independent), 



+ 



(2^ 



Moreover, the momentum independent contribution appears to diverge in the — 
limit. Indeed, 

1 f°° 

Vra (momentum independent) = / dril + e"^™""^ + (1 + 74)'Tir)x 

2m^ JO 







' A A^^ 
1 + ^ + 7^ 



2 9 92 

+ 



A9A 9A2 



D(r,A)+ 



1 r°° 

Jo rf^(l + e-''"" + (l-74)mr) X 



(29) 



A 



x/ c/A<^ 1--1 



2A 2A3 



3r 



+ 



r 



4^ 



A 

r 



1 



A A2 A3 

+ 



2 9 92 



+ 



D{t,X), 



2r 2r2 4r3 / \^A9A SA^ 
and for large r (r ^ T^) one has the following asymptotics for expression (^91) : 



(a+(l + 74) + a-(l-74))^, mT, > 1 



(momentum independent) 



/9 er r 
m i„ 



mT„ < 1, 



(30) 



where a+, a_ and /5 are some coeffitients of order unity depending on the explicit form of 
function D{t,X). It is clear from expression (^) and asymptotics (^) that for the case 
mTg ^ 1 Vp is indeed a corection to the leading regime (p3|), (^41). For mT^ <^ 1 this 
correction is larger than the leading interaction, and we conclude in such a way that the 
embarrasing regime (|2^) is not potential, and, moreover, for such light quarks one should 
turn to the full system of equations (|ll]) and ([12|) . 
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The term proportional to the binding energy is simply 



and does not cause many additional problems, it is small comparing to the confinement 
term if eoTg ^ 1, as it was already found in 

The effective Schroedinger equation is obtained from Dirac equation ([T9|) by the stan- 
dard Foldy-Wounthuysen (FW) reduction, and displays a lot of pleasant features for 
mTg ^ 1. For the upper Dirac component of the quark wave function the resulting 
Hamiltonian is 

Hfw = ni + ^ + eE{r) + VLs{r) + £Af(r) + Vsi{r), (32) 
where SE{f) is the standard confining interaction 

eE{r) = 2 dr dXD^r, X){r - X), (33) 



and 

VLsir) = --^ r dr f dXD{T, X) f 1 - - V (34) 
2m"^r Jo JO V / 

Expressions (^) and (34) were obtained from equation (p!9|) in |^, and coincide with 
the ones given by the Vacuum Background Correlators Method [3-6]. As it was shown in 
P], form (|34D for the spin-orbit force is equivalent to the form obtained by Gromes [§] in 
terms of Wilson loop expectations. 

As for r ^ Tg 

Ee = err 

and 

Vls = , 
4m2r ' 

these expressions mimic scalar confinement (0) at large distances. Nevertheless, confining 
interaction (|2^) is not a scalar one, at large distances it is 5/6 scalar and 1/6 vector, 
and is due to the electric correlator ^^44 only as one should expect for the nonrelativistic 
particle. 

There are three sources for the spin-orbit force (^4]). One piece comes from the FW 
reduction of confinement (^3]), it is purely magnetic and gives 



v!:S = -^rdrrdXDir,X) l-^] - 
^■^ 3mV Jo Jo \ r3 / r»Tg Qmh 
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The second one stems from the FW reduction of ^yk term and contains both electric and 
magnetic contributions yielding 



/1^ (T r°° f 




r>T9 12m?r 

The third piece is from k'^ term. It is purely magnetic and its share is 



3m^r Jo Jo \ 2r ) < 

Magnetic confinement 



oo rr I ') \ \^ 



JO 



£M(r) = 2r / rfre-^'"^^ / dXD^r, A) - - - + — , (35) 



3 r 3r3 



behaves as at large distances and it is suppressed for large mTg. 

There are, of course, spin-independent corrections of the standard Darwin term type, 
which are O (^y^:^ at large distances. Another correction is given by electric momentum 

dependent part of Vp , which is O , and the binding energy correction ( PT| ) contributes 
at the same order of magnitude. Nevertheless, the main spin-independent correction 
comes from the electric part of (momentum independent): 

Vsj(r) = - I dTT I dX 



m Jo Jo 1^7-2 



which behaves as at large distances, competing with magnetic confinement (|35|) . 

The results of the suggested approach reproduce well-known fomulae for leading con- 
finement (|33| ) and spin-orbit ( p4|) forces. Our derivation does not appeal to the Wilson 
loop approach and effective QCD string at large distances, but, as model (0) for the 
QCD vacuum is compatible with the area law, the salient features are the same. A kind 
of a string is developed connecting quark and antiquark, and this string is the minimal 
string of the Vacuum Background Correlators Method [4-6] , or the flux tube with gluonic 
degrees of freedom in the ground state M, lll|, as far as we neglect the A^^^ contributions 



to the field strength correlator (plsl). Nevertheless, as we deal with full Dirac Hamiltonian 
([TP|), where both nonlocality in time direction and Zitterbewegung are included, we are 
able to disclose new important corrections (|35D and (^). 

We have demonstrated that, apart from naive nonrelativistic condition m ^ -y/a and 
more sofisticated condition eoTg -C 1 0, another requirement, mTg ^ 1, is needed for 
the potential-type description of heavy quark dynamics to be valid. To what extent 
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this requirement is indeed new? In the Nc — ^ oc hmit correlators ([T5|) are given by 
the pure Yang-Mills theory, with single nonperturbative mass scale. On the other hand, 
model ( ]T5| ) contains two dimensional parameters, correlation length Tg and -D(O), which is 
proportional to the gluonic condensate and is related to the string tension a via equation 
(|T6|). So it is not surprising that two phenomenological quantities, ^/o■ and T~\ should 
be of the same order of magnitude, and indeed they are. The commonly accepted value 
of the string tension gives ~ QAGeV , while lattice measurements |jT2[ give for the 
correlation length value Tg^ ~ IGeV, so the requirement mTg ^ 1 does not bring drastic 
changes into our understanding of heavy quark dynamics. Much more interesting are new 
corrections (|35|) and (0), which are potentially large and even at large distances depend 
on the explicit form of vacuum correlation function D. 

We acknowledge useful discussions with Yu.A.Simonov. 
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